Competition between singlet and triplet pairings in Na x . C0O2 • y H 2 
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We discuss the pairing symmetry of a cobaltate superconductor Na^CoCb • y H2O by adopting an 
effective single band model that takes into account the e' g hole pockets, as discussed in our previous 
paper [to appear in Phys. Rev. Lett.] Here we consider the off-site repulsions in addition to the 
on-site repulsion considered in our previous study. We show that the spin-triplet /-wave pairing 
proposed in our previous study is robust to some extent even in the presence of off-site repulsions. 
However, /-wave pairing gives way to singlet pairings for sufficiently large values of off-site repulsions. 
Among the singlet pairings, i-wave and extended s-wave pairings are good candidates which do not 
break time reversal symmetry below T c in agreement with the experiments. 



I. INTRODUCTION 



Discovery of superconductivity in a hydrated sodium 
cobaltate Na K Co02 -yH^Oi has attracted much attention 
recently. This material is of special interest since it is a 
3d-electron system having a quasi-two-dimensional lat- 
tice structure, which resembles the high T c cuprates. On 
the other hand, it differs from the cuprates in that the Co 
ions form a triangular lattice, and also there are several 
bands intersecting the Fermi level, all of which are away 
from half-filling. Another interesting point concerning 
this material is that some experiments suggest possible 
occurrence of an unconventional superconductivity^*^ 
In particular, there has been a debate as to whether the 
pairing occurs in the spin singlet channel or in the triplet 
channel. An unchanged Knight shift across T c found 
in some experiments suggests spin-triplet pairing, 3 while 
the Knight shift does decrease below T c in other experi- 
ments at low magnetic fieldi^ As for the occurrence of 
broken time reversal symmetry, several [iSR experiments 
show that such a possibility is rather small&S 

There has also been a singlet-triplet debate theoret- 
ically. Some theories propose singlet pairing ^iii^ii&ii, 
others triplet pairing mechanisms ^i^i^l^i^ Naively, 
presence of possible ferromagnetic spin fluctuations 
reported in some experiments^ may support spin- 
fluctuation-mediated triplet pairing scenario. Couple of 
years ago, however, two of the present authors showed 
that T c of spin triplet superconductivity due to ferromag- 
netic spin fluctuations, if any, should in general be very 
lowi^i A similar conclusion has been reached by Mon- 
thoux and Lonzarich. 22 This is because the triplet pairing 
interaction mediated by spin fluctuations is proportional 
to x/2 (in contrast to 3%/2 for the singlet case), where 
X is the spin susceptibility, while the effective interaction 
that enters the normal self energy is proportional to 3x/2, 
so that the suppression of pairing due to the normal self- 



energy overpowers the pairing effect. We then showed 
that this difficulty for spin-fluctuation-mediated triplet 
pairing may be eased in systems having 'disconnected 
Fermi surfaces', where the nodal lines of the gap can 
run in between the Fermi surfaces to open up a full gap 
on the Fermi surfaces^ Since pairing with high angu- 
lar momentum is degraded by the presence of gap nodes 
on the Fermi surfaces, opening a full gap may result in 
an enhanced p a i r i n g,2i 2 4i2£i 2 iii2L 2 £ As an example, we 
considered the Hubbard model on a triangular lattice, 
where the Fermi surface becomes disconnected into two 
pieces centered around the K point and the K' point (see 
Fig^a)) when there are a small number of holes (or elec- 
trons, depending on the sign of the hopping integral). Us- 
ing fluctuation exchange (FLEX) approximation 29 ' 30,31 
and solving the linearized Eliashberg equation, we have 
shown that a finite T c can be obtained for spin-triplet 
/-wave pairing, where the nodes of the gap run be- 
tween the Fermi surfaces m& We have also confirmed this 
conclusion 3 ^ using dynamical cluster approximation, 33 
which is a non-perturbative approach. 

In a recent study^ we have looked at the band calcula- 
tion results of Na.i;Co02 34 from this point of view, where 
we find that the pocket-like Fermi surfaces (hole pock- 
ets) near the K and K' points (as in the bottom figure of 
Fig^b)), originating from the band having e' g character 
to some extent in the notation of refl34l (denoted by the 
thick line in Fig^a)), are disconnected in a similar sense 
as in the triangular lattice, namely the nodes of the /- 
wave gap do not intersect the Fermi surfaces (although 
they do intersect the large Fermi surface around the T 
point originating from the bands having ai g character). 

In fact, a strong motivation to focus on the e' g hole 
pockets is the presence of van Hove singularity (vHS) 
near the K point. Namely, as can be seen from the band 
structure shown in Fig|2 there exist saddle points at 
points denoted as SP, where the density of states (DOS) 
takes a large value. Since this large DOS lies close to the 



Fermi levels, it is likely that the band structure around 
the K and K' points strongly affects the low energy prop- 
erties of this material. In particular, ferromagnetic spin 
fluctuations may arise due to this high DOS near the 
Fermi level. 



In refH we have adopted a model where we separate 
out the portion of the bands which has e' g character to 
some extent, and in particular focus only on the upper e' g 
band, which contributes to the formation of the pocket 
Fermi surfaces and to the large DOS at the vHS, while 
neglecting the lower e' g band, which has only small DOS 
near the Fermi level due to the linear dispersion at the 
band top. Although this portion of the band, shown by 
the thick curves in Fig|^a), is disconnected between M 
and r points due to a\ g -e' g hybridization, it originally 
comes from a single band having e' character, as can 
be understood from the inset of FigJ^a), where a tight 
binding band dispersion is given for a case when large 
aig-e' g level offset is introduced. We have found that the 
thick portion of the band in Fig^a) (apart from the 
missing part between M and T points) can be roughly 
reproduced by a single band tightbinding model on a tri- 
angular lattice with hopping integrals up to fourth near- 
est neighbors. Namely, the dispersion of the tight bind- 
ing model on an isotropic triangular lattice with only 
nearest neighbor hoppings takes its maximum at the K 
point (FigQJa)), while the band top moves towards the 
r point when the fourth neighbor hopping is introduced 
(Fig^b)), resulting in a pocket like Fermi surface that 
lies between the K and T points when small amount of 
holes are present. Moreover, the band structure around 
the K point resembles that of the actual material (com- 
pare FigQJb) and the thick curve in Fig|2J), so the vHS 
due to the saddle points SP in Figgis also reproduced 
by this effective model. As for electron-electron interac- 
tions, we considered a model with only the on-site re- 
pulsion, where we found using FLEX that spin triplet 
/-wave pairing indeed dominates and has a finite T c i& 
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FIG. 1: Band dispersion (top) along with the Fermi surfaces 
(schematic) for a small number of holes and the /-wave gap 
(bottom) are shown for a tight binding model on a triangular 
lattice with (a) only nearest neighbor hopping, and (b) when 
fourth nearest neighbor hopping = 0.2 is introduced. + — 
denote the sign of the gap function, while the dashed lines 
represent the nodal lines. 



II. FORMULATION 
A. Model 

In standard notations, the Hamiltonian considered in 
the present study is given in the form, 

H= Y1 (/'■•'-''"' + h - c -) + u n ^ ni i 

ij , <t i 

+Vi ^2 n icr n j<j' +V2 ni a n ja > (1) 

(ij)',a;a' (ij)",cr,a' 



In the present study, we continue to study the sin- 
gle band model having hoppings up to fourth nearest 
neighbors, but this time with repulsive interactions up to 
second nearest neighbors to investigate the effect of the 
off-site repulsions. We find that (i) spin triplet /-wave 
pairing is robust to some extent even in the presence of 
off-site repulsions, but (ii) for sufficiently large off-site re- 
pulsions, singlet pairings can dominate over triplet pair- 
ings. 



where Uj — t\, t-i, t^, and are the nearest, second 
nearest, third nearest and fourth nearest neighbor hop- 
ping integrals, respectively^ and U, Vi, V2, represent on- 
site, nearest neighbor, second nearest neighbor electron- 
electron repulsions, respectively. The summations Yl(ijy 
and y)ijjy, are taken over pairs of nearest and next near- 
est neighbors, respectively. We define the hole density as 
rt/j = 2 — n, where n is the electron band filling (number 
of electrons/number of sites), rih is restricted to a small 
number in the present study because the e' g hole pockets 
are small. We take —t\ — 1 as the unit of the energy 
throughout the study. 
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FIG. 2: (a) A schematic plot of the band structure of 
Na^CoCh following the calculation results of refl34l. The thick 
line denotes the portion of the upper part of the bands which 
has e' g character to some extent. Bottom inset: dispersion of 
a three band tight binding model, where large ai g -e' g level off- 
set as well as third nearest neighbor hoppings between same 
orbitals is introduced in addition to the hoppings considered 
m refpia. (b) Brillouin zone is shown in the extended zone 
scheme. The points denoted as SP are saddle points since the 
energy of the thick band at SP is lower than those around 
the gray circles, while it is higher than those around the open 
circles. 



B. Method 

To make the calculation feasible at low temperatures 
even in the presence of off-site repulsions up to second 
nearest neighbors, here we adopt a random phase approx- 
imation (RPA) approach. The singlet (V s ) and triplet 
(V ) pairing interactions are given within RPA as fol- 
lows. 



V s {q) = U + V{q) + -U 2 X s(q) 



--(U + 2V(q)) 2 Xc (q) 



-\{U + 2V(q)fxM) 



(2) 



(3) 



Here, V(q) is the Fourier transform of the off-site repul- 
sions. Xs and Xc are the spin and charge susceptibilities, 
respectively, which are given as 



xo(q) 



xM) = 



1 - U X o(q) 

xo(q) 



l + (U + 2V(q)) Xo (q)' 



(4) 



Here xo is the bare susceptibility given by 

f( £ P+q) ~ f( £ P ) 



xo(q) 



with Ek being the energy dispersion measured from the 
chemical potential and f(e p ) is the Fermi distribution 
function. We take 128x128 fc-point meshes in the actual 
numerical calculations. 

To obtain the onset of the superconducting state, we 
solve the gap equation within the weak-coupling theory, 



AA(fc) = -^V<*(fc 



tanh(/? £k ,/2) 

2£fc' 



The transition temperature Tc is determined by the con- 
dition, A = 1. In the present approach, lu dependence as 
well as the self energy correction is neglected. Although 
this approximation is quantitatively insufficient, we be- 
lieve that the present approach suffices for the purpose 
of discussing the competition (relative tendency toward 
pairing) among various pairing symmetries. 



C. Possible pairing symmetries 

Pairing symmetries are classified according to the irre- 
ducible representations of the point group Dg, as shown 
in Fig|3 For the cases we considered, the /-wave pair- 
ing having B2 symmetry turned out to have small values 
of A, so in the following sections we consider s, p, d, /, 
j-wave pairings having A\, E%, E2, B±, A2 symmetry, 
respectively. 




s (A,) d (E 2 ) i (A 2 ) 




p (EO /(BO /(B 2 ) 



FIG. 3: Possible pairing symmetries of the present model 
classified according to the irreducible representations of D6. 
The dashed lines represent the basic nodal lines. 
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III. RESULTS 

A. Case without off-site repulsions 

We first show results for a case without off-site repul- 
sions, where we take ti = £3 = 0, £4 = 0.18, i%h = 0.18, 
U = 2.4 and Vi = V 2 = 0. As seen in Fig[j] the /- 
wave pairing strongly dominates over other pairing sym- 
metries, which is consistent with our previous study^ 
The reason for this is: (i) the spin fluctuations are nearly 
ferromagnetic (see Fig|SJa)), and (ii) the nodes of the f- 
wave gap do not intersect the Fermi surfaces as seen in 
FigEl 




0,6 0.024 0.028 0.032 0.036 0.04 



T 

FIG. 4: The eigenvalue of the gap equation plotted as a 
function of temperature for each pairing symmetry, — tz = 
0, U = 0.18, U = 2.4, Vi = V 2 = 0, n h = 0.18. 




FIG. 5: Contour plot of the gap functions for each pairing 
symmetry. The dashed lines represent the nodes of the gap, 
and the white rings are the hole pockets (only two out of the 
six are shown for clarity). Parameter values are the same as 
in FigUl and T = 0.03. All the gap functions are plotted at 
T = 0.03 throughout the study. 



B. Case with moderate nearest neighbor repulsion 

We now turn on the off-site repulsions. As a case where 
a moderate nearest neighbor repulsion exists, we take 
V\ = 0.7 and V2 = 0.4, while other parameter values are 
taken to be the same as in section UlI Al We see in FigEl 
that the /-wave pairing still dominates over the others, 
but it is suppressed compared to the case without off-site 
repulsions. On the other hand, A of singlet i,d 7 and s- 
wave pairings, along with the singlet gap functions shown 
in Fig[7| are barely affected by the off-site repulsions. To 
understand the reason why /-wave is suppressed while 
singlet pairings are unaffected by off-site repulsions, we 
show in FiglSJthe pairing interaction in momentum space. 



Here we define the "spin part" and the "charge part" of 
the pairing interaction as U 2 Xs(q) and V(q) — ^(U + 
2^(q)) 2 Xc(<7), respectively. The spin part peaks near 
the r point, while the charge part takes a broad negative 
minimum value around the M point. Thus, the effect of 
the pair scattering between two pockets that is mediated 
mainly by the charge part cancels out when there are 
nodes on the pockets (Fig^Ja)), while the sign change 
of the /-wave gap is unfavorable in the presence of the 
charge part, which has a negative sign (FigEJb)). 



C. Case with large nearest neighbor repulsion 

For a larger nearest neighbor repulsion, where we take 
the parameters values to be the same as in section UlI Bl 




T 

FIG. 6: A plot similar to FigHjwith Vi = 0.7, V 2 = 0.4, and 
other parameter values taken to be the same. 



except that V\ — 1.0, we find that /-wave pairing gives 
way to singlet pairings, as shown in Fig llOl Among the 
singlet pairings, which are found to closely compete in a 
wide range of parameter values, i-wave pairing surpris- 
ingly has the largest T c for the present set of parameter 
values. Note that /-wave is no longer full-gaped on the 
Fermi surfaces due to the additional nodes as seen in 
FigHU 



IV. DISCUSSION 
A. Intuitive picture in real space 

The reason why triplet pairings are suppressed by 
off-site repulsions while singlet pairings are barely af- 
fected can intuitively be understood in a real space pic- 
ture. As seen in Fig|S] and Fig|7| the triplet pairing 
gap functions have relatively small number of nodes, 
which shows that these pairing takes place at short dis- 
tances in real space. In particular, /-wave gap in the 
absence of off-site repulsions (Fig|SJl essentially has the 
form sin(^Y-k x + ^ky)— sin(fc y )+sin(— ^k x + ^k y ), which 
can be rewritten, apart from a constant, as 

exp(ik ■ a) — cxp(ik ■ b) + cxp(ik • c) 
— exp(— ik ■ a) + exp(—ik ■ b) — cxp(— ik ■ c), (6) 

where a = (^,|), b = (0,1), c = (-^,|), and k = 

{k X t ky ) • 

Since ±a, ±b, and ±c are the coordinates (in units 
of the lattice constant) of the six nearest neighbors, the 
/-wave pairing in the absence of off-site repulsion is es- 
sentially a nearest neighbor odd-parity pairing. By con- 
trast, the singlet pairings have many additional nodal 
structures compared to the basic nodal structure shown 
in FigEJ This means that the pairs are mainly formed 
at large distances in real space. In fact, the result that A 
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FIG. 7: Gap functions for the parameter values of Fig|S| 



is unaffected by V\ and V2 is consistent with the picture 
that the singlet pairs are formed at distances farther than 
second nearest neighbors. 

The reason why singlet pairs are formed at large dis- 
tances while triplet pairs are formed at short distances 
can intuitively be understood from a nearly ferromag- 
netic spin alignment in real space as shown in Fig ll2l 
Namely, electrons at close distances have nearly parallel 
spins, while the electrons can have antiparallel spins only 
when they are separated far away. 

In order to reinforce this real space intuitive picture, 
we look into a case where the spin susceptibility peaks 
at positions more closer to the T point, namely when 
the spin structure is more purely ferromagnetic than in 
the cases considered in sections IlII Bl and IIII CI We take 
i4 = 0.15 and ri/j = 0.1, which moves the vHS and the 
Fermi level towards the band edge, thus making the spin 
structure more ferromagnetic as seen in Fig ll3f a). 

As shown in Fig lMf al. the triplet pairings dominate 
strongly over singlet pairings, and among the triplet pair- 




FIG. 8: Contour plots of the spin (a) and the charge (b) part 
of the pairing interactions for the parameter values of FiglUl 
The spin part takes its maximum value ~ 25 at a short wave 
but finite vector Q s , while the charge part takes its minimum 
value ~ —3 at the M point. 



ings, p-wave has the largest T c . This result can again be 
understood from the above real space picture. Namely, 
since the spin susceptibility peaks around the T point, 
the wave length of the spin alignment is large, so that 
the singlet pairing distance is large from Fig^| In fact, 
the number of nodes in the singlet pairing gaps are found 
to increase compared to those in sectior HIl CI as seen in 
FigEH This makes the singlet pairing less competitive 
against triplet pairings. Moreover, since the wave length 
of the spin alignment is large, not only electrons residing 
at nearest neighbors but also those separated to some ex- 
tent have parallel spins, so that spin triplet pairing can 
take place at relatively large distances, thereby circum- 
venting the effect of off-site repulsions. It can be seen 
from Fig^l that the /-wave gap has additional nodes 
close to the Fermi surfaces suggesting that the pairing 
takes place not only at nearest neighbors but also at 
larger distances. Pairing at large distances, and thus the 
presence of nodes near the Fermi surfaces, degrades /- 
wave pairing, so that p-wave pairing dominates. 

Even if we decrease V\ down to V\ — 0.7, p-wave still 
dominates over /-wave as seen Fig lMf b 1 ). in contrast to 
the case discussed in sections IIII Bl This is because in 
the presence of purely ferromagnetic spin fluctuations, 
/-wave pairing takes place at large distances even for 
relatively small off-site repulsions, so that the additional 
nodes are located close to the Fermi surfaces. This can in 
fact be seen by comparing Fig ll5l and Fig llGI where the 
additional nodes in the /-wave gap do not move away 
so much from the Fermi surfaces even if we reduce Vy 
from 1.0 to 0.7. This is in contrast with the case of Fig[7| 
and Fig |12l where the additional nodes in the /-wave 
gap clearly moves away from the Fermi surfaces (i.e., the 
nearest neighbor pairing component becomes large) as 
we reduce V\. 
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(a) 




FIG. 9: The reason why singlet pairings (here we take the 
i-wave gap in Fig0 as an example) are barely affected by 
off-site repulsions (a), while /-wave pairing is suppressed (b). 
Q c is the wave vector at which the charge part of the pairing 
interaction takes its maximum absolute value. 



B. Peculiar Disconnectivity of the Fermi surfaces 

The reason why i-wave pairing, which has a very large 
angular momentum, may dominate in the presence of off- 
site repulsions can be found in the disconnectivity of the 
Fermi surface peculiar to the present system. Namely, 
the Fermi surfaces in the present case are disconnected 
in a way that they do not intersect the T-M lines nor 
the hexagonal Brillouin zone edge, both of which are the 
nodal lines of the i-wave gap. Thus, the number of nodal 



FIG. 10: A plot similar to Fig© with Vi = 1.0 and other 
parameter values taken to be the same. 



lines that intersect the Fermi surfaces is similar among 
the three singlet pairing symmetries, extended s, d, and 
i-waves, as seen in Fig llll thereby making the competi- 
tion subtle. In fact, we have looked into various sets of 
parameter values and found that the three singlet sym- 
metries tend to have close values of A. 

To show more clearly that the subtle competition be- 
tween i-wave and other singlet pairings is due to the dis- 
connectivity of the Fermi surfaces peculiar to the present 
system, we compare the above results with those for cases 
with only nearest neighbor hopping, namely with a sim- 
ple triangular lattice. In a simple triangular lattice with 
rih < 0.5, the Fermi surfaces are disconnected in a man- 
ner that they do not intersect the T-M lines, but they 
do intersect the Brillouin zone edge. The absence of t± 
moves the van Hove singularity away from the band top, 
so in order to have a moderate T c as in the cases above, 
we either have to take larger values of n,/,, or increase the 
values of the electron-electron interactions. Here we take 
i%h = 0.4 and the interactions the same as in sectior llll CI 
in Fig ll7f ah while in Fig llTf b') we take the same but 
take the interactions 1.5 times larger than in sectior llll CI 
It can be seen that in both cases, A of the i-wave pair- 
ing is much smaller than that of the d- and the s-wave 
pairings. This is because the Fermi surfaces intersect the 
Brillouin zone edge, i.e., the nodes of the i-wave gap, as 
seen in FiglTland FiglBI 



C. Comparison with the case with 
antiferromagnetic spin fluctuations 





FIG. 11: Gap functions for the parameter values of Fie llOl 



singlet 



triplet 
< 



%/Q 



FIG. 12: An intuitive real space picture of singlet and 
triplet pairings when nearly ferromagnetic spin fluctuations 
are present. 



In the present study, we have looked into a model 
where nearly ferromagnetic spin fluctuations and charge 
fluctuations due to off-site repulsions coexist. It is inter- 
esting to compare the present situation with the case 
where antiferromagnetic spin fluctuations and charge 
fluctuations coexist. We have previously proposed that 
spin-triplet pairing is possible in a quarter-filled, quasi- 



1D organic material (TMTSF)2X^S*21 where supercon- 
ductivity lies right next to spin density wave (SDW) in 
the pressure-temperature phase diagram, so that the spin 
fluctuations should be antiferromagnetic. Naively, an- 
tiferromagnetic spin fluctuations would favor spin sin- 
glet d-wave pairing, while experiments suggest triplet 
pairingi22ii£ Our proposal is that if 2fcjr (= tt/2 because 
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FIG. 13: The spin (a) and charge (b) part of the pairing 
interaction for t% = t s = 0, t 4 = -0.15, U = 2.4, Vi = 1.0, 
P2 = 0.4, n n = 0.1. The spin part takes its maximum value 
~ 35 at a vector Q a ~ 0, while the charge part takes its 
minimum value ~ — 3 at the M point. 



the band is quarter-filled) charge fluctuations coexist 
with 2kp spin fluctuations, as suggested from an exper- 
imental fact that CDW actually coexist with SDW 41 i 42 , 
spin triplet '/-wave' like pairing, which is essentially a 
fourth neighbor pairing having a gap form of sin(4fc x ) 
(this is called /-wave in the sense that the gap changes 

sign as H 1 1 — along the Fermi surface), may become 

competitive against singlet 'd-wave' like pairing, which is 
a second neighbor pairing with a gap form cos(2k x ). In 
momentum space, this close competition arises because 
(i) the number of nodes of the /-wave gap on the Fermi 
surface is the same with that of the d-wave because of 
the disconnectivity of the Fermi surfaces due to quasi- 
one-dimensionality, and (ii) the pairing interactions due 
to spin and charge fluctuations (apart from the first or- 
der terms) have the form V s = §P sp — ^V c h for singlet 
pairing and V 1 = — \V sp — j^ch for triplet pairing, where 
V sp and V c h are contributions from spin and charge fluc- 
tuations, respectively, so that the absolute values of V s 
and Vt becomes comparable when V sp ~ P c h- A similar 
theory has been proposed by Fuseya et alm& Recently, 
we have confirmed this scenario microscopically by ap- 
plying RPA to an extended Hubbard model with off-site 
repulsions up to third nearest neighbors^ Here, a key 
parameter is the second nearest neighbor repulsion, due 
to which the 2kp charge fluctuations arises. 

In real space, this f-d competition can be understood 
as follows. When there are purely antiferromagnetic spin 
fluctuations, second neighbor singlet pairing is likely to 
occur, as seen from Fig l2(JI If we now turn on the second 
nearest neighbor repulsion, the second neighbor singlet 
pairing is degraded, and a fourth neighbor triplet pairing 
becomes competitive. Since this intuitive picture in real 
space is rather general, the tendency that spin singlet 
pairing, which is favored by antiferromagnetic spin fluc- 
tuations, gives way to triplet pairings in the presence of 
sufficiently large off-site repulsions, should hold in vari- 
ous situations. In fact, such a singlet-triplet crossover due 
to off-site repulsions is seen also in a half-filled quasi-lD 



(a) 




0.028 0.032 0.036 0.04 



T 



(b) 




0.028 0.032 0.036 0.04 



T 

FIG. 14: Plots similar to Fig© for ti = t 3 = 0, U = 0.15, 
U = 2.4, V 2 = 0.4, n n = 0.1. Pi = 1.0 (a) and Pi = 0.7(b). 



system^ a model for Sr 2 Ru04 on a 2D square lattice^ 
and a model for Na^CoC^ on a 2D triangular latticed 

Now, it is interesting to point out that the above ten- 
dency is exactly the opposite of what has happened in the 
present case with nearly ferromagnetic spin fluctuations 
: triplet pairings formed at relatively short distances give 
way to singlet pairings formed at farther distances when 
sufficiently large off-site repulsions are introduced. Thus, 
the effect of the off-site repulsion totally differs depend- 
ing on whether the spin fluctuations are ferromagnetic 
or antiferromagnetic. It is also important to note that in 
both cases, the disconnectivity of the Fermi surfaces can 
assist the pairings at large distances (singlet in the ferro- 
magnetic case and triplet in the antiferromagnetic case), 
which have larger number of nodes in the gap function. 
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FIG. 15: Gap functions for the parameter values of Fig lMf a). 



D. Validity of focusing on the pockets 

In the present study, we have focused only on the e' g 
Fermi pockets. Although it is necessary to consider a 
multiband modeUS in order to strictly take into account 
the effect of the Fermi pockets, our single band approach 
is supported by a recent multiband FLEX study^ where 
a similar /-wave pairing is found to dominate when spin 
fluctuations are not too purely ferromagnetic. 

Experimentally, the Fermi pockets are not observed 
in ARPES experiments^^ up to date. However, since 
these experiments are done for materials with relatively 
large Na content, (i.e., large number of electrons in CoO 
layers) it is likely that the Fermi level lies above the 
e' g bands4& In fact, an experimental result that maxi- 
mum T c is reached only when the content of Na decreases 
sufficiently 4 — can be considered as an indirect support for 
scenario in which the e' g band plays an important role for 
superconductivity. 




FIG. 16: Gap functions for the parameter values of Fig ll4f b'l . 



E. Pairing symmetry in the cobaltate 

Let us finally discuss the pairing symmetry in the 
cobaltate in light of the present results along with the 
existing experiments. We have seen that "full-gapped" 
/-wave proposed in ref[3 is robust to some extent 
even in the presence of off-site repulsions. Existence of 
nodes in the gap function is suggested in a number of 
experiments^^ which may seem to contradict with the 
/-wave scenario. In the actual Na^OoC^, however, there 
is also the large Fermi surface around the T point, which 
is not taken into account in the present study. Then we 
can propose a possible scenario where the 'active' Fermi 
surfaces for superconductivity are the pockets, while the 
large Fermi surface is a 'passive' one, in which supercon- 
ductivity is induced by the 'active' band via interband 
interactions. Since the nodes of the /-wave gap intersects 
the large Fermi surface, this can account for the experi- 
mental results suggesting the existence of gap nodes. In 
fact, it has been shown that a /iSR data can be explained 




by taking this view£ 

On the other hand, singlet pairings are found to domi- 
nate over triplet pairings in the present model when spin 
fluctuations are not purely ferromagnetic and the off-site 
repulsions are sufficiently large. Among the singlet sym- 
metries, z-wave is another good candidate that competes 
with /-wave since it does not break time reversal symme- 
try below T c . We have shown that i-wave may dominate 
over other symmetries, but the parameter range where 
this symmetry dominates is not so wide as that for the 
/-wave. Although extended s-wave is also consistent with 
almost all the existing experiments, we could not find a 
parameter regime where this symmetry dominates. How- 
ever, considering the fact that the competition within the 
singlet pairings is very subtle, and that the present model 
only roughly takes into account the electronic structure 
of the actual material, we believe the possibility of ex- 
tended s-wave pairing cannot be ruled out. 

Since d-wave and p-wave pairings are likely to result 



in a broken time reversal symmetry and/or a full gapped 
state (unless an accidental situation occurs) below T c due 
to the two-fold degeneracy, we believe that the possibility 
of these pairing states are less likely. We may conversely 
say from the present results that the spin fluctuation in 
the actual cobaltate is not purely ferromagnetic and/or 
the off-site repulsions are not too large, so that p-wave 
pairing does not dominate. 

At present, it is not clear experimentally whether the 
pairing occurs in the singlet or in the triplet channel be- 
cause Knight shift experiments are done on polycrystal 
samples. In order to settle this singlet-triplet debate, 
Knight shift experiment for high quality samples, or al- 
ternatively, tunneling spectroscopy studies based on a 
newly developed theory for triplet superconductors pro- 
posed by one of the present authors^ is required. Thus, 
our viewpoint is that possibility of neither singlet nor 
triplet pairings cannot be ruled out experimentally at the 
present stage. Then, the bottom line obtained from the 
present study, combined with the existing experimental 
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facts, is that triplet /-wave pairing is most likely, al- 
though there remains a possibility of singlet i-wave pair- 
ing, and also a possibility of extended s-wave pairing 
cannot be ruled out. This ambiguity is a direct con- 
sequence of the very point of the present study. Namely, 
the present system is highly peculiar in that anisotropic 
pairings having large angular momentum, such as /- and 
i- waves, are good candidates due to the disconnectiv- 
ity of the Fermi surfaces. The fact that these "unusual" 
pairing symmetries are good candidates makes the pair- 
ing competition subtle. 



FIG. 19: Gap functions for the parameter values of Fig ll7f bi. 



SUMMARY 



To summarize, we have studied an effective single band 
model for the hydrated sodium cobaltate focusing on the 
e' g hole pockets, where we take into account off-site re- 
pulsions up to second nearest neighbors. We find that 
/-wave pairing proposed in our previous study is robust 
to some extent in the presence of off-site repulsions. For 
large off-site repulsions, /-wave gives way to singlet pair- 
ings that does not break time reversal symmetry below 



triplet 



ffl_J* tit* 



singlet 
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FIG. 20: An intuitive real space picture of singlet and triplet 
pairings when antiferromagnetic spin fluctuations are present. 
Note that when Q = 2k F = tt/2 as in the case of (TMTSF) 2 X, 
arrows in this figure are separated by two lattice spacings. 
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